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ROUNDED FRONTS AND BLUNT BASES 
By Stanley H. Scher and James S. Bowman, Jr. 

SIMMARY 


An Investigation has Been made of the stability of several Bodies 
of revolution having fineness ratios smaller than 1.0 and having rounded 
fronts and Blunt Bases . Free-drop tests were made in the vertical 
Langley 20-foot free-spinning tunnel to determine the Behavior of the 
various Bodies in descent. The results indicated that the stability of 
the Bodies as determined from observed motions decreased as the Body 
shapes varied from a hemisphere to a near-paraboloid. Free-oscillation 
tests were ma de for some of the Bodies and the static stability param- 
eter Cn^ was found to Be Indicative of the relative Behaviors noted 

during the free-drop tests. In addition to the determination of sta- 
bility from the free -Body tests, a method for calculating the approxi- 
mate stability of the subject type of Body was developed. Results 
obtained By this calculation method agreed well with the experimental 
results . 


INTRODUCTION 


Recently, interest has Been expressed in the Behavior in descent of 
various Bodies of revolution with rounded fronts and Blunt Bases and with 
fineness ratios less than 1.0. Inasmuch as the theories, such as those 
presented in references 1 and 2, which are available for calculating the 
stability of Bodies of revolution are not applicable to Bodies of such 
small fineness ratio, experimental free-drop tests to determine the inher 
ent Behavior of several such Bodies have Been made. The tests were made 
in the vertical air stream of the Langley 20-foot free-spinning tunnel, 
and the Body shapes included a hemisphere, a near -paraboloid, and four 
intermediate Bodies . Some of the Bodies (those which were stable) , 
while free in the air stream, were oscillated in a single plane and 
their degrees of stability were determined from observed damping times 
and periods of the oscillations . In addition, a method for calculating 
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the approximate static stability of the subject types of bodies is 
developed and applied to the bodies . 


SYMBOLS 


x 

r 

L 

lY 


a 

Cl Vo< 

s 

q 

p 

T l/2 


axial distance from bow of body to any body station 
local body radius at any station x 
body length 

moment of inertia about lateral body axis through center of 
gravity perpendicular to longitudinal body axis 

free-stream velocity 

angle of attack (body-axis inclination relative to free-stream 
flow direction) 

drag coefficient of body at an angle of attack of 0° ; Weight 

qb 

area of base of body 
free-stream dynamic pressure 
period of oscillation 

time for oscillation to damp to one -half applitude 


dC 


m 


^ “ da 


K 


su mm ation of aerodynamic derivatives contributing to moment 
caused by pitching velocity of body 

Mass 


K-i 


relative density factor , 

p Volume 

nondimens ional radius of gyration corresponding to Iy, k y /L 
differential operator, 

nondimens ional time parameter based on length, Vt/L 
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k Y 


radius of gyration corresponding to Iy 


t time 


P 


mass density of air 


9 

d F 
dM 




M 


angle "between tangent to body contour and longitudinal body 
axis in plane of symmetry in which a is measured 

force on element (positive when upward along transverse axis 
or rearward along longitudinal axis) 

pitching moment about center of gravity due to force on element 
(positive when sense of moment is nose up) 

axial distance from bow of body to center of gravity 

moment-arm distance (in plane perpendicular to longitudinal 
body axis) from the axis to center of gravity of projection 
of surface area of half the hollow cylinder element (used 

-r, the value for a semicircumference, representing half of 

jt 7 

an infinitely thin hollow cylinder )j see fig. 2 

pitching-moment coefficient , — — ^ 

* q. Volume 

pitching moment about center of gravity due to integration of 
forces on all body elements (positive when sense of moment 
is nose-up) 


c d<x,=0 o local longitudinal-force coefficient at an angle of attack 
of 0° at any x-station based on body circumference 

c do=90° local transverse-force coefficient at an angle of attack of 90 ° 
at any x-station based on body diameter 

circular -cylinder-section drag coefficient based on cylinder 
c diameter 

H ratio of the drag coefficient of circular cylinder of finite 

length to that for cylinder of infinite length 



NACA RM L52L08 


k 


Subscripts : 

t transverse direction (perpendicular to longitudinal body axis) 

l longitudinal direction (parallel to longitudinal body axis) 

av average 


MODEL AND APPARATUS 
Model 


Profile sketches and over-all dimensions of the bodies of revolution 
used in the investigation are presented in figure 1. The ordinates of 
the bodies axe presented in table I. As can be seen in figure 1, the 
bodies included a hemisphere, a near -paraboloid, and four intermediate 
bodies . The six bodies shall sometimes be referred to in the remainder 
of this paper as bodies A through F, by which letters they are identified 
in figure 1 and table I . Each body of revolution was symmetrical about 
its longitudinal axis . For the experimental part of the investigation, 
the bodies used were made homogeneously of balsa segments and had lac- 
quered surfaces. Mass and design information for the bodies are presented 
in table II, the center-of-gravity locations as listed in the table being 
the natural ones for the respective homogeneous bodies. 


Wind Tunnel 

The Langley 20-foot free-spinning tunnel, which has an airspeed 
range up to about 95 feet per second, was used, for the wind-tunnel test 
phase of the investigation. This tunnel has a vertically rising air 
stream and was convenient to use in investigating experimentally the 
behavior in descent of the bodies of revolution. 


EXPERIMENTAL TESTS 
Free-Drop Test Methods 


In the free-drop tests, each body was held with its longitudinal 
axis in an attitude somewhere between vertical and an angle of attack 
of about 20° and was released to float in the tunnel air stream. Visual 

observations were made of the behavior of the body in the tunnel air 

stream. Descent velocity was recorded as the airspeed necessary to hold 

the body at a given level in the tunnel. For those bodies which descended 
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in a nose-down vertical attitude, the drag coefficient of the hod y at an 
angle of attack of 0° was obtained by dividing the weight of the body by 
the free-stream dynamic pressure at descent velocity and by the projected 
area of the body in a plane perpendicular to the longitudinal body axis . 
The results were analyzed to determine the effects of varying the body 
shape on the behavior of the bodies in descent. Reynolds numbers of 
the free-drop tests, based on lengths of the longitudinal axes and on 

the rates of descent of the bodies, ranged from about 2 .5 x 10** to about 

7.5 x 10** . Reynolds numbers based on cross-body components of velocity 
(V sin a) and on body diameters at various stations along the longitudinal 

axis were for the most part between 1.0 x 10** and 3 *0 x 10**, although 
the lower limit naturally approached zero for forward body stations with 
very small diameters or when the angle of attack of a body was very 
small. 


Free-Drop Test Results 

The results of the free-drop tests are presented in table III. 

After being released to float in the tunnel air stream, the hemispherical 
and near -hemispherical bodies of revolution, bodies A and B in figure 1, 
quickly assumed a nose-down vertical attitude with little or no hori- 
zontal travel in the tunnel. Bodies C and D also trimmed generally in 
a nose-down vertical attitude, although for body C there were occasional 
oscillations of the order of ±5° from the vertical, approximately about 
a lateral axis through the center of gravity, and for body D, there was 
some slight horizontal travel in the tunnel. Body E trimmed in a nose- 
down vertical attitude part of the time but also tended to make erratic 
falling-leaf motions which eventually evolved to a continuous end-over- 
end tumbling motion approximately about a lateral body axis through the 
center of gravity while traveling horizontally in the tunnel. Body F, 
a near-paraboloid, descended with a tumbling motion similar to the one 
just described. The results Indicated that the stability in descent of 
the homogeneous balsa bodies A to F decreased as the body shapes became 
less like that of a hemisphere. Although this decrease in stability may 
have been due in part to more rearward (percentagewise) center-of -gravity 
locations, a comparison of the stability of bodies C and F, which had 
identical center-of -gravity locations, indicates a definite effect of 
shape on stability. A brief additional test was made in which lead 
ballast was attached to the front of body F, causing the center of gravity 
to move forward far enough so that the body trimmed in a nose -down 
vertical attitude . If sufficient turning moment was applied at launching 
to start any of the bodies rotating, the body would descend in the tunnel 
with a tumbling motion similar to that described for body E. The causes 
and geometry of these tumbling motions are not considered in this paper, 
inasmuch as such a study is beyond the scope of the present Investigation. 
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Free-Oscillation Test Methods 

As mentioned earlier, it was noticed during the free-drop tests that 
the motion of four of the subject bodies of revolution (A, B, C, and F 
with its center of gravity moved forward) as they floated in the air 
stream of the wind tunnel appeared to be substantially a single -degree - 
of -freedom motion approximately about a lateral axis through the center 
of gravity of the body. Based on the assumption that such a motion was 
present, motion pictures were made as the bodies were forced to oscillate 
by means of an externally applied moment, and the period of the oscilla- 
tions and the time to damp to one-half amplitude were obtained from the 
film records. By using the values of P and T values of C^, as 

well as of K (representing the summation of the derivatives contributing 
to the moment caused by the pitching velocity of the body) , were deter- 
mined as shown in the appendix. 


Free-Oscillation Test Results 

The values of K and obtained for the four bodies from 

oscillation during the free drop Eire listed in table IV. As illustrated 
for body B in the appendix, the term containing K in the root of the 
equation of motion when applied to each of the bodies was found to be so 
small that it had only a negligible effect on the values of (static 

stability) obtained. The degrees of static stability of the various 
bodies as determined from the free -oscillation tests showed qualitative 
agreement with the relative degrees of stability indicated by the motions 
observed during the free-drop tests. 


STABILITY CALCULATIONS 
Method of Calculating Stability 


As air flows ar ound a body of revolution of the subject type when 
inclined at an angle of attack, the velocities over the bottom surface 
are different from those on the top half, and these velocity differences 
cause unbalanced pressures and forces which may cause a pitching moment 
on the body. An attempt is made here to develop a convenient method for 
estimating the pitching moments of the bodies without the necessity of 
making a detailed study of pressures and forces at various points on the 
front or base surfaces. Consideration is given, by a strip-analysis 
method, to forces on the body normal to the surface profile in the plane 
of symmetry in which a is measured, the forces being taken as functions 
of the free-stream velocity components in that normal direction. In using 
this approach, cognizance was taken of some two-dimensional results pre- 
sented in reference 3 which indicate that the force on a circular cylinder 
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in the defined normal direction is a function of the component of free- 
stream velocity in that direction. The (frictional) flow parallel to 
the surface profile is disregarded. 

It is convenient to consider a body element of length dx (see 
fig. 2) as being the wedge-shaped front of a hollow right -circular- 
cylinder element with a length of L - x, a thickness of dr, and pro- 
jected surface areas in longitudinal and transverse (parallel and 
perpendicular, respectively, to longitudinal body axis) body planes 
of 2r dx and 2«r dr. When the body of revolution is at an angle of 
attack, the components of free-stream velocity normal to the surface at 
the bottom and top points of the body element dx in the plane in which 
a is measured may be expressed, respectively, as V sin (0 + a) and 
V sin (0 - a) . On the base of the body (rear surface of the hollow 
cylinder) directly behind the bottom and top points on the front, the 
free-stream velocity components normal to the surface may each be 
expressed as V cos a. For convenience, let us consider the unbalanced 
force on the hollow-cylinder element due to these different normal 
velocity components in terms of separate transverse and longitudinal forces. 
At the bottom and top points, respectively, of the front of the element, 
the transverse force components would be proportional to 


— V 2 sin 2 (0 + a) cos 0 
2 


and 


£• V 2 sin 2 (0 - a) cos 0 
2 


and the longitudinal force components would be proportional to 


— V 2 sin 2 (0 + a) sin 0 
2 


and 

£• V 2 sin 2 (0 - a) sin 0 
2 

(vectors shown in fig. 2) . On the rear surface of the hollow cylinder 
there would be no transverse force components and the longitudinal force 
components at both the bottom and top points of the element would be pro- 
portional to £■ V 2 cos 2 a. The pitching moment about the center of gravity 
of the body due to the transverse and longitudinal forces on the element is 
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dM = dFt(xm - x) + cLFjr e 


( 1 ) 


Expanding equation (l) by use of the foregoing relationships, and including 
the local force coefficients cd a= 9Q° 821(1 c d a _oo Sives 

dM = c^ a _^ 0 ° jq sin 2 (0 + a) cos 0 - q sin 2 (0 - a) cos ©] (x m - x) 2r dx - 

c d a=0 c t sin 2 (0 + a) sin 0 - q cos 2 aj2r e rtr dr + 

c da=0° & sin 2 (0 - a) sin 0 - q cos 2 aJ 2r e nr dx 

or 

dM = 2qc da _ 90 o cos 0 |sin^(0 + a) - sin 2 (0 - a)J (x m - x)r dx ■ 

2 JtqCd a _ 0 o sin 0 |sin 2 (0 + a) - sin 2 (0 - a)Jr e r dr 

It can be shown by expanding and collecting terms that 

|sin 2 (0 + a) - sin 2 (0 - a)j = sin 2a sin 20 
and, since r e , 


dM = 2q sin 2a c. _ cos 0 sin 20 (x m - x) r dx - 

<^= 90 ° 

4q sin 2a c d a _QO sln ® sln 2 ® ( 2 ) 

The moment due to the forces over the vhole body may be obtained by 
integrating the incremental moments over the length of the body. 

In the present investigation, the value of c D a= oo (table III) 

determined from the free-drop tests of each body was used as c^.qo 

for all body elements of that particular body. In regard to values for 
c ■, it was believed that the present purpose of calculating the 

a txc =90 
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approximate stability of the bodies could be achieved by using the value 
of the product ^^(cos 0) a v as c d a= 90 ° for 8,11 elements of a 

particular body. The quantity qc^ was originally intended (ref. 4) 

for use in determining the average section drag coefficient of a circular 
cylinder of finite length and the (cos 0) av term is intended as an 

attempt to compensate for the fact that T)c ^ does not include the 

effects of body curvature along the length of the present bodies . 

Although values of tjc^ were obtained (ref. 4) only at a cross-flow 

Reynolds number of 8.8 x lcA, they are assumed herein to be applicable 
over the Reynolds number range of the present investigation, in which 
range c^ c has a constant value of 1.2. A value of q was obtained 

for each body by extrapolating the curve shown in figure 3 which is based 
on the results of the aforementioned experiments on circular cylinders. 

A value of (cos 0) av was obtained for each body by plotting a curve of 

cos 9 against body station, integrating the area under the curve, and 
then dividing by the body length. 


In the present investigation, integrations of the products of the 
variables in the two terms of equation (2) were performed graphically 

2 sin 2 <xc^_^qO 

and multiplied by the respective factors - 


Volume 


and 


k sin 2 ac 


for an angle-of-attack range from 0° to 90° to obtain the 


Volume 

pitching-moment coefficient over this range, 


Results of Stability Calculations 

The results of the calculations are presented in figure k and in 
the tabulated values of in table V. As would be expected from the 

nature of the equation used in calculating the pitching moment, the 
curves for C m against a are symmetrical about a maximum value at an 
angle of attack of 45° and return to a value of zero at an angle of attack 
of 90°. The results indicate that all the bodies except body F are 
statically stable, and the results are therefore in qualitative agreement 
with the wind-tunnel test results discussed earlier in this paper. For 
the various statically stable bodies, the slopes of the curves show a 
general agreement between the relative degrees of stability as obtained 
experimentally and as calculated; that is, bodies A and B are most stable, 
body C less stable, and so forth. Good quantitative agreement was 
obtained between the values from the calculation method and those 

obtained from the oscillation techniques, as may be seen by comparing 
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tables IV and V. It Is of interest that the results of the calculations 
for body F with its center of gravity moved forward to 50 -percent body 
length indicate stability; this result is in agreement with the experi- 
mental results. Table V includes a breakdown of the calculated resultant 
stability of each body in terms of the contributions to caused by 

the respective transverse and longitudinal forces. As can be seen from 
the table, the contributions due to the longitudinal force for each bo y 
was a stabilizing moment, whereas the contribution due to the transverse 
force was destabilizing except for body F with its center of gravity 
moved forward. 

It is of interest that, if an average value of tased on 

the six measured values obtained for the bodies of the present investi- 
gation had been used in calculating the stability of each body, the 
changes in would be within about ±20 percent of the values calcu- 

lated. Thus, it appears that, if the method developed for calculating 
the stability of the subject bodies is applied to other bodies of the 
same general type and no measured value of ca^o is available, an 

average or estimated value of c^^o can be used and should allow 
qualitative determination of the stability of the body with a reasonable 
degree of reliability . 


CONCLUDING REMARKS 


An investigation has been made of the stability of several bodies 
of revolution having fineness ratios smaller than 1.0 and having rounded 
fronts and blunt bases . Free-drop tests were made in the vertical Langley 
20-foot free-spinning tunnel to determine the behavior of the various 
bodies in descent. The results indicated that the stability of the bodies 
as determined from observed motions decreased as the body shapes varied 
from a hemisphere to a near -paraboloid. Free -oscillation tests were made 
for some of the bodies and the static stability parameter C^ was found 


to be indicative of the relative behaviors noted during the free-drop 
tests. In addition to the determination of stability from the free-body 
tests, a method for calculating the approximate stability of the subject 
type of body was developed. Results obtained by this calculation method 
agreed well with the experimental results. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Iangley Field, Va. 
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APPENDIX 

METHOD OF DETERMINING STABILITY OF BODY BY OSCILLATION TECHNIQUE 


The applicable nondimens ional single -degree -of -freedom equation of 
motion for the bodies of revolution is 

2PKY 2 D 2 a - & Da - C^a = 0 


from which the characteristic equation is 


2UK y 2 D 2 - | D - C^ = 0 


or 


D 2 . K D _ J25* 0 

4WCy 2pKy 


The roots are 


The formulas 


D = ± 

8pK Y 



+ 


Cn bc 

2pK Y 2 



and 


p _ 2* L 

Imaginary part of root V 

were used to determine K and C^. order to illustrate the use of 

this technique, work done for body B is shown. From an extrapolation of 
the plot in figure 5, which was prepared from motion-picture records of 
the oscillations of body B, it was determined that 

T-^ = 1.17 sec 
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Then, into the formula 

T = 1.17 sec 


0.693 (8nK y 2 ) l 
K V 


the following constant values for "body B axe substituted 


Therefore, 


L = 0.1475 ft 
V = 1+6.5 ft/sec 

p = 11+6.8 

k y 2 = 0.1953 

K = -0.1+305 


where the negative sign is inserted with the quantitative value of K 
"because the damping opposes the oscillatory motion. 

From figure 5, it was also determined that 

P = 0 .35 sec 


Then, from the formula 



and "by using the constants and K for body B it was found that 


Cn^ = -0.186 


where the minus sign is inserted with the quantitative value of 
because the static stability opposes the oscillatory motion. 
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TABLE IV.- VALUES OF K AMD OBTAINED FROM OSCILLATION TECHNIQUE 


Body 

K 

c ma 

per radian 

A 

1.242 

- 0.272 

B 

.430 

-.186 

C 

1.21 

-.128 

ap 

1.38 

-.149 



Center of gravity of "body F moved forward to 
50-percent body length. 
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Figure 3 .- Ratio of the drag coefficient of a circular cylinder of finite 
length to that of a cylinder of infinite length as a function of the 

fineness ratio. Data taken from reference A. Reynolds number, 8.8 x 10 
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Figure 5.- Variation of amplitude of oscillation with time for body B 
as obtained from motion-picture records . 





